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IN THE TITLE: 

Please change the title to read: 

METHOD AND DEVICE FOR USING A MULTI-CHANNEL MEASUREMENT 
SIGNAL IN SOURCE MODELLING 
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IN THE SPECIFICATION: 
Page 2, 1 st full paragraph 

Conventional minimum norm estimates involve inherent problems such as slowness 
of calculation and susceptibility to noise. For example, in the case of an L2 norm, one needs 
an inverse matrix of matrix G, whose element (i, j) contains the inner product of the 
switching lead fields of the ith and jth measurement sensor, so one must calculate these inner 
products for each pair of sensors. The switching lead field is so determined that the signal 
measured by a sensor is the projection of the current distribution for the switching lead field 
of the sensor in question. The noise problems are due to the fact that matrix G calculated for 
the sensors is susceptible to noise, so in the calculation of its inverse matrix, regularisation is 
needed in the practical situations. 

Page 2, 2 nd full paragraph on page. 

Regularisation methods, such as the br e ak off truncation regularisation of the singular 
value decomposition, usually are non-intuitive, and usually also to be solved for each case 
specifically. A regularisation of the wrong type may lead to an erroneous modelling result. 

Page 3, last paragraph on page continuing to top of page 4 

The basic idea of the invention is that because the computation of the inner products 
of the sensor fields is hard and difficult using a conventional set of sensors, it is worth using a 
special set of sensors, whose switching lead fields are orthogonal and, if possible, to be 
analytically computed. In principle, this can be implemented by a suitable physical set of 
sensors. As a suitable physical set of sensors is, however, often quite difficult to manufacture, 
it is, in most of the cases, advantageous to use virtual sensors computationally generated from 
a conventional set of sensors, i.e. the measurement signals are converted into other ones by a 
suitable conversion so that they correspond to the signals that the virtual measurement device 
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would have measured. At the same time, it is possible, if necessary, to eliminate the signals 
associated with external interference. This conversion has been described e.g. in patent 
application FI20030392, which is incorporated herein by reference. After the conversion, the 
source modelling is performed in an optimal manner using the basis vector components of the 
signal space instead of the actual measurement signals. One substantial feature of the 
invention is that after the conversion, the source model need not be any more regularised. 
Page 4, last paragraph on page. 

According to the invention, a multi-channel measurement signal corresponding to 
each measurement sensor is converted into the signals of a predetermined set of virtual 
sensors, and the current distribution of the object being examined is determined by depth r 
from the signals of the set of virtual sensors in a beforehand selected orthonormal function 
basis. In that case, the estimation of a current distribution is fast and robust. Further, to 
achieve the set of signals corresponding to the set of virtual sensors, a multi-pole 
development is calculated from a multi-channel measurement signal. A multi-pole 
development can be calculated in two ways: by taking into account the magnetic fields 
emitted by sources outside the object being measured, or by ignoring them. 

Page 5, 1 st paragraph on page. 

Advantageously, as the orthonormal function basis, a basis with the following form is 
selected: 

jCr) = 2 5> ta /-,(r)i/m(<9,0>), 

/=0 m=-/ 

wherein f-(r) fi(r) is a radial function to be freely selected andXi m (0,<p) Xi m (0, <p) is 
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so-called vector spherical harmonic. In that case, it is possible to place the function basis into 
a current distribution equation, and the coefficients of the current distribution are analytically 
solved based on the equation: 



c lm -r> M ; 



_0 



1-1 

c lm = r, M lm 



]r' +2 Mr)dr 

_0 



wherein y { is a constant associated with order 1 and R is the radius of the 

sphere to be examined. Advantageously, function f(r) yj(r) is used to adjust the depth 

weighing of a current distribution model. 

Page 5, last paragraph on page continuing to top of page 6. 

According to the invention, the processing means include a conversion module for 
converting a multi-channel measurement signal corresponding to each measurement sensor 
into the signal of a predetermined set of virtual sensors; and calculation means for 
determining the current distribution of an object being examined or for calculating by depth r 
from the signals of the set of virtual sensors in a beforehand selected orthonormal function 
basis. In one embodiment, the calculation means are arranged to calculate a multi-pole 
d e velopment expansion from a multi-channel measurement signal. 

Page 7 5 3 rd full paragraph on page. 

In the following section, the mathematical background and grounds of the invention 
are described. When the magnetic fields are converted into coefficients M lm = a lm + ib lm 

associated with the basic solution r~ (M) Ylm(0,<p) of the Laplace equation, wherein i is an 
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p p -> 

imaginary unit, they can be expressed by means of the current distribution KrT A r ) in 
spherical coordinates (r, 0, p) , whereby they are of the form: 

M lm -y { \r l hm{8,cp).j{r)dv ^ 

V 

M lm =r,jr l Xlm{9, <p) • J(7)dv QX 

v 

wherein the integration is performed over the volume being examined, y t is a constant 
p -> 

associated with order 1 and Xlm\d,<p) Xlm{0, (p) is so-called vector spherical harmonic. 
Entire page 8 continuing to top of page 9. 

This form is apparent can be derived e. g. from publication "Multipole expansions of 
electromagnetic fields using Debye potentials", C. G. Gray, American Journal of Physics, 
Vol. 46, pp. 169-179, 1978. The expression mentioned above is of the switching lead field 
form, wherein the switching field of the multi-pole coefficient M lm is of the form: 

Li m {r) + r x Xi m (9,<p) (3± 

Lt m (7) + r l Xi m (0,<p) £21 

On the other hand, the vector spherical harmonics form by depth r an orthonormal 
basis, so with the depth in question, the current distribution can be presented in the function 
basis in question: 

J (r) = g g e lm f x (r) Xlm(0 9 <p), 

J(r) = lL £c lm f,(r)Xlm(0,<p), £31 

1=0 m=~l 

wherein f t (r) is some radical function. 

10195538.1 



Docket No. 032700-8 
Serial No. 10/572,570 
Page 7 



When as the volume to be examined, spherical volume is selected, by placing 



the previous expression into the equation (1), the coefficients of the current distribution can 
be analytically solved: 



wherein y x is a constant associated with order 1 and R is the radius of the sphere to be 
examined. The previous equation (4) indicates that the coefficients of a current distribution 
model presented in an orthonormal basis can be solved based on coefficients M lm in a 
completely trivial manner using analytical expressions without any kind of regularisation. 
This is computationally very fast and numerically stabile. Function f(r) f t (r) ls freely 

selectable, and can be used to adjust the depth weighing of a current distribution model. 
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